Heavy-quark contribution to the proton's magnetic moment 
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We study the contribution to the proton's magnetic moment from a heavy quark sea in quantum 
chromodynamics. The heavy quark is integrated out perturbatively to obtain an effective dimension- 
6 magnetic moment operator composed of three gluon fields. The leading contribution to the matrix 
element in the proton comes from a quadratically divergent term associated with a light-quark tensor 
operator. With an approximate knowledge of the proton's tensor charge, we conclude that a heavy 
sea-quark contribution to the proton's magnetic moment is positive in the asymptotic limit. We 
comment on the implication of this result for the physical strange quark. 
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1. Introduction. Naively, the proton is composed 
of three valence quarks bound by the strong interac- 
tion. However, quantum chromodynamics (QCD), the 
theory of strong interactions, implies a much more com- 
plex structure: The proton also contains a sea of virtual 
gluons and quark- antiquark pairs. These virtual com- 
ponents of the proton have measurable consequences on 
the proton's macroscopic properties, such as its mass and 
spin. In this article, we study the contribution to the 
magnetic moment of the proton from a heavy sea quark 
with mass ttlq 3> Aqcd, the strong interaction scale. 

Our study is motivated by the recent results from a se- 
ries of remarkable experiments carried out at MIT-Bates, 
Jefferson Lab, and Mainz 1]. They all seem to indicate, 
within experimental errors, that the strange quark sea 
makes a positive contribution to the proton's magnetic 
moment in the standard convention that leaves out the 
quark's electric charge factor. On the other hand, most 
models, along with some lattice-based calculations, pre- 
dict a negative strange quark contribution to the proton's 
magnetic moment |2| . Therefore, there seem to be a con- 
tradiction between experimental data and theoretical re- 
sults. A resolution of this contradiction might help us to 
better understand the influence of the strange quark on 
the macroscopic properties of the proton. 

A direct calculation of the strange quark contribution 
to the proton's magnetic moment is a challenging task, 
mainly because the strange quark is neither too light nor 
too heavy compared to Aqcd: Its mass lies just at the 
border between these two regimes. This is unfortunate 
because various successful theoretical methods have been 
developed for these two limits only. On the other hand, it 
might be useful to treat the strange quark mass as a vari- 
able, and to study the limiting cases where the strange 
quark is either very light or very heavy. 

There is strong evidence that a light sea-quark con- 
tributes negatively to the proton's magnetic moment, al- 
though a rigorous proof does not yet exist. In this arti- 
cle, we focus on the opposite limit: We consider a heavy 
strange quark and evaluate its contribution to the pro- 



ton's magnetic moment. By integrating out the heavy 
quark perturbatively, we obtain an effective magnetic 
moment operator made of pure gluon fields. We find that 
the leading contribution of the matrix element is related 
to the nucleon's tensor charge and is positive. We thus 
conclude that a sea-quark contribution to the proton's 
magnetic moment depends non-trivially on the sea-quark 
mass: It must change sign at some critical value of the 
sea-quark mass. We argue that this critical value could 
be smaller than the physical mass of the strange quark. 

2. Effective gluonic operator from heavy quark cur- 
rent. Consider a proton of momentum p and polarized in 
the z-direction, |pt). The- z-component of its magnetic 
moment is defined as 



( P 1\yd 3 r(rx j em )>t) 

(pTIpT) 



(1) 



where j cm is the electromagnetic current. In order to 
carry out the spatial integration, it is useful to consider 
the off-diagonal matrix element and to take the forward 
limit after integrating over r 



Vv = -\ (v> x (p'tl J(0) H)) ; 



q=0 



(2) 



where q = p' — p is the momentum transfer. 

Now consider the contribution to the magnetic momen- 
tum from a heavy quark with electromagnetic current 
3q = Ql^Q- (The quark also contributes to the mag- 
netic moment through wave functions. However, this is 
conventionally regarded as a part of the light-quark con- 
tribution.) To lowest order in 1/raq, its contribution to 
the form factor is given by the diagram shown in Fig. 1. 
When ttiq Aqcd , the quark loop can be integrated out 
perturbatively. The heavy-quark electromagnetic current 
is matched to a leading effective operator consisting of 
pure gluon fields [3| 



j^Q = C{m Q ) d a T^ a (m Q ) + 



(3) 
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where ellipses denote higher order contributions in I/tuq, 
and C(ttiq) = g 3 (thq) / ((4tt) 2 Abrnq) to leading order in 
the strong coupling g(rriQ), and 

Tfj, a = 14Tr G^{G a \G Ta }-5T r G aT {G aT , G afl } , (4) 

where G M1/ is the strong field-strength tensor. Notice that 
the tensor T MQ is antisymmetric and that the effective 
gluon operator is renormalized at the scale tuq. 




FIG. 1: Lowest order heavy sea quark contribution to the 
proton's magnetic moment. 

The matrix element of the heavy-quark current can 
thus be written as 

(p'Tb'MQbT) = C(m Q )(p']\d a T m \p]) 

= iC{m Q )q a {p]\T^\p]) . (5) 

From Eq. (2), the heavy-quark contribution to the pro- 
ton's magnetic moment is then given by 



6$ = C(m Q ) (ptl^xbt) , (6) 



with 



abc 



7{E a -B b ) E cz -2(E a -E b -B a -B b ) B c 



(7) 

where d abc are the SU(3) C coefHcients, and E a and B a 
are the chromo-electric and chromo-magnetic fields, re- 
spectively. We have followed the convention in the liter- 
ature which leaves out the electric charge, eQ, factor. As 
is obvious from Fig. 1, the effective operator is exactly 
the same as that obtained for light-by-light scattering at 
lowest order in quantum electrodynamics (QED) 

3. Top quark contribution to the magnetic moment of 
the Ab baryon. To understand the nature of the gluonic 
matrix element in the proton, let us consider a similar 
case in QED. At order a^ m , the electron's anomalous 
magnetic moment receives a contribution from the muon 
through light-by- light scattering. In the limit m^/m* ^> 
1, this quantity has been analytically calculated in |5j: 



from fi 



I «cm \ J / m e 
V 7T / V m„. 



-C(3) - — 
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This is a striking result: The leading power dependence 
on the muon mass is rather than 1/ro^ as implied 

by a naive heavy-m^ expansion similar to (4-8). The 
discrepancy is resolved if one realizes that the effective 
operator TJ X (obtained by replacing the gluon fields by 
photon fields in Eq. (5)) is renormalized at scale m^. 
Power counting indicates that the matrix element of T3 X 
in the electron state is quadratically divergent, and the 
divergence is cut-off naturally by the scale to m . Since the 
matrix element of the effective operator has dimension 3, 
it must behave as ((e|e) = (27r) 3 <5 3 (0)) 



{e\ T yx\ e ) ~ m l m e > 



(9) 



where the m e dependence comes from the chirality flip 
of the electron. This demonstrates that the matrix ele- 
ment of the effective operator is dominated by quantum 
fluctuations at the scale of the heavy fermion mass. 

The QED result can be directly applied to an interest- 
ing case in QCD. Consider a spin-1/2 Ah baryon made of 
a heavy bottom quark and two light up and down quarks 
coupled to zero spin and isospin. To leading order in 
1/rrifc, the spin of Aj, is carried by the heavy 6-quark. The 
magnetic moment of A b from the top quark sea can be 
calculated just like in the QED case. Including a color 
factor of Cd — ^Z a b c d abc /4& = 5/18, the leading order 
contribution in Aqcd/w?,, mb/m tl and a s is 



s tA„ = A* A t (^) C d ( 



m t 



:C(3) 
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(10) 



where fiA b — eth/2rriA b c. We thus conclude that the top 
quark contributes positively to the Af, magnetic moment. 

4. Gluonic matrix element in the proton. The above 
example is very instructive for the behavior of the effec- 
tive gluonic operator matrix element in the proton: Its 
dependence on the renormalization scale tuq and Aqcd 
must exhibit the form, 



{P]\T V X\P]) 



am Q A QC D 



(11) 



where the possible logarithmic dependence on itiq is ne- 
glected. The a-term depends quadratically on tuq and 
is therefore dominant. It receives contributions mainly 
from quantum fluctuations at scale tuq. Power count- 
ing of the diagrams involved in i|ll|) shows that this term 
comes from single quark contributions. Diagrams involv- 
ing two or three quarks from the proton cannot produce a 
term quadratic in toq. The 6-term is entirely determined 
by physics at the non-perturbative scale Aqcd- 

Wc first focus on the a-term since it dominates the ma- 
trix element (|llfl • Because it involves contributions from 
both large momentum flow and scale Aqcd , we factorize 
them by matching the dimension-6 operator T^ v to a set 
of dimcnsion-4 ones 



'i 1 



(8) 



C, 



(12) 
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where 9% must be gauge-invariant second-order antisym- 
metric tensors with odd charge parity, and ellipses denote 
higher dimensional operators. Since we are interested in 
the forward matrix element, total derivative operators 
are ignored. Such gauge invariant dimension-4 operator 
cannot be constructed with pure gluonic fields. There 
are, however, two such operators made of quark fields, 

e^ a ^{iD alp - iDp 7a )iP, m4,o^i> . (13) 

Using the QCD equation of motion (i 1/) — m)ip — 0, the 
first operator can be reduced to the second one. There- 
fore, we find that 



(14) 



where a = i(a,j/7r) 3 Cd[3/2£(3) — 19/16] comes from the 
matching calculation for a single quark state, and the 
sum runs over the light-quark flavors. 

It is now straightforward to calculate the leading con- 
tribution to the magnetic moment, 



Sfip = C (mQ)m,Qa(m u Su + maSd) 



(15) 



where Su and Sd are the tensor charges of the proton de- 
fined as (PSi^ia^^lPS) = 2Sip{S"P u - S U P»), with 
P and S the proton's four-momentum and spin, respec- 
tively . The above equation is one of the main results 
of this paper. Using the PDG values m u = 2.8± 1.3 MeV 
and rrid — 6 ± 2 MeV Q and the tensor charge from the 
MIT bag model Su = 1.2 and Sd = -0.3 @, we find that 
S/i® is positive. Lattice calculations at higher renormal- 
ization scales yield a similar ratio Su/Sd 8] which does 
not change the conclusion. 

The strange quark also contributes to Eq. (16) through 
the proton wave function, but its contribution is sup- 
pressed compared to that of the light up and down 
quarks. The scenario we study here is that the strange 
quark is the fictitious heavy quark Q. The corresponding 
mqSQ term in Eq. (16) can be calculated by matching 
Qa^ v Q to the gluon operator T^ v and then evaluating 
the leading contribution from the matrix element T MI/ in 
the proton state. The result is that although rriQdQ is 
a leading power, it is suppressed by a^mq) relative to 
the up and down quark contributions and is therefore 
significantly smaller. 

5. An example of non-perturbative contribution. The 
subleading non-perturbative contribution is hard to cal- 
culate without a specific model of the proton. Here we 
present quark model calculations where the excitations of 
the quarks and the non-linear gluonic interactions are ne- 
glected. In this Abelian approximation, the gluon fields 
consist of eight copies of the electromagnetic field obey- 
ing the linear Maxwell equations. The chromo-electric 
fields are given by E a (f) = ¥-E(r*) = -^V$(r), with 



d 3 n 



pin) 

\r- f!| 



(16) 



where the index a refers to its color, p(r) = (a\p(f)\(3) 
denotes the charge density matrix element between two 
states, \a) and |/3), and A a are the Gell-Mann matri- 
ces. Similarly, the chromo-magnetic fields are given by 
B a (r) = ^-B(r) = ^-V x A(r) with 



d 3 ri 



J in) 
r- fi| 



(17) 



where j(r) = (a\j(f)\{3) is the quark current density ma- 
trix element. 

We compute the contribution to the magnetic moment 
due to a single quark polarized in the +z-direction, and is 
in the ground state of a quark model potential. The elec- 
tric and magnetic fields can be calculated using the wave 
functions that describe the quark fields in the ground 
state of the model. We use the non-relativistic (NR) 
quark model and the MIT bag model 0, . The domi- 
nant contribution to the proton magnetic moment is 



5(1® — C(mg) / dr r 2, m(r) 



(18) 



m{r) 



7E^ • B-[-[ E- 



TT G TT 



2E^ • E\x B 



TT -°TT 



2 ( -Btt • B\-[ B 



TT -°TT 



-Bn ■ P> B- 



HT ^TT 



; ReB n ■ B n B z u 



The numerical evaluation of these integrals results in 
Sfi^/n N C(m Q ) ~ 2.2 fm- 4 and 1.3 fin" 4 in the NR 
quark model and in the MIT bag model, respectively. 
If naively applying this result for the physical strange 
quark, with m s = 110 MeV and a s (m s ) ~ 1, we get 
Sfip ~ 0.1e s ^jv, where /in = eh/2Mc is the nuclear mag- 
neton. 

In both models, the Sfj,® is therefore positive. However, 
this result is reached for different reasons. In the NR 
model, the B field is generated by the color dipole of 
the quark and hence is correlated with its spin. The 
triple- B term can be neglected because the quark is non- 
relativistic. When averaging over the direction of the B 
and E inside the proton, the term 7(E ■ B)E Z dominates 
over 2E 2 B Z . In the MIT bag, on the other hand, the color 
current is generated from the orbital motion of the quark. 
The direction of motion is correlated with the spin. Here 
the triple-S term is non-negligible. After canceling the 
negative E 2 B Z term, the result is again positive. 

It would be nice, of course, to also calculate contribu- 
tions from intermediate excited states, particularly the 
contribution with excitation energy of order toq. Be- 
cause of the intermediate state sum, the single-quark ma- 
trix clement dominates the 6-term. Such a calculation is 
beyond the scope of this paper. In the NR quark model, 
it can be shown that all excited state contribution is pos- 
itive. For the MIT bag, on the other hand, we have only 
verified this for the first few excited states. 
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6. Mass-dependence of the sea quark contribution to 
the proton's magnetic moment. Thus it appears that a 
heavy sea quark contributes positively to the proton's 
magnetic moment. On the other hand, there are indica- 
tions, that a light sea quark contributes negatively. 

One sensible picture for the light quark sea is the me- 
son cloud model 0,^3- ^ the strange quark is light, the 
kaon is then a Goldstone boson. The proton has a certain 
probability to dissociate into a kaon, which contains an 
anti-strange quark, plus a A, which contains a strange 
quark. The kaon carries a unit of orbital angular mo- 
mentum and its contribution to the magnetic moment is 
negative. The A's spin and magnetic moment are dom- 
inated by the strange quark and the latter is positive. 
However, in the virtual dissociation, the A is polarized 
mainly against the proton spin. Thus the strange quark 
contribution to the magnetic moment is also negative. 
Therefore the total contribution from light strange and 
anti-strange quarks is negative. 

There is also an indirect lattice QCD calculation for 
the nucleon sea contribution to the magnetic moment of 
the proton |l2^ . 

^" 1o °p = Sfi s p - 0.3(1) , (19) 

where (5/ip~ loop is the up (or down) quark "sea contri- 
bution" defined in the same sense as the strange quark 
sea. Thus unless the real strange quark contribution is 
larger than 0.3fi n (the recent world data average gives 
<5/4(Q 2 - O.lGeV 2 ) = (0.28 ± 0.20)e s ^ N 13]), the light- 
sea contribution shall be negative. 




FIG. 2: A possible picture for a sea-quark contribution to 
the proton's magnetic moment as a function of the sea-quark 
mass. 

The ineluctable consequence of these results obtained 
in the light and heavy quark limits is that the sea-quark 
contribution to the proton's magnetic moment vanishes 
at least once at some critical value of the sea-quark mass. 
How large is this critical mass? A plausible answer is that 
its value is slightly above the masses of the up and down 
quarks, i.e. around 10 to 20 MeV. Indeed, this is where 
the chiral behavior seems to rapidly set in, as can be seen 
from chiral extrapolations of lattice data 12]. If this is 
correct, the strange quark contribution to the magnetic 
moment of the proton is positive. 



7. Conclusion. In this paper, we advocate the study of 
the contribution to the magnetic moment of the proton 
from a heavy-quark sea. We show that the leading contri- 
bution can be calculated in perturbative QCD combined 
with the quark tensor charge and that the contribution is 
positive. If one believes that the contribution is negative 
when the quark mass is small, the magnetic moment has 
a nontrivial dependence on the sea-quark mass. There- 
fore the sign of the strange quark contribution to the 
proton's magnetic moment is very sensitive to the QCD 
dynamics. 
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